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Abstract: We explore the possibility of a new dark matter candidate in the supersym- 
metric type III seesaw mechanism where a neutral scalar component of the Y= triplet 
can be the lightest supersymmetric particle. Its thermal abundance can be in the right 
range if non-standard cosmology such as kination domination is assumed. The enhanced 
cross-section of the dark matter annihilation to W + W~ can leave detectable astrophysical 
and cosmological signals whose current observational data puts a lower bound on the dark 
matter mass. The model predicts the existence of a charged scalar almost degenerate with 
the dark matter scalar and its lifetime lies between 5.5 cm and 6.3 m. It provides a novel 
opportunity of the dark mater mass measurement by identifying slowly-moving and highly- 
ionizing tracks in the LHC experiments. If the ordinary lightest supersymmetric particle 
is the usual Bino, its decay leads to clean signatures of same-sign dilepton and di-charged- 
scalar associated with observable displaced vertices which are essentially background-free 
and can be fully reconstructed. 



Contents 



[l]. Introduction 

[2]. Triplet spectrum in supersymmetric type III seesaw 



|3|. Relic density and indirect detection limits 



LHC signatures 



|5|. Conclusion 



1. Introduction 

Dark matter (DM) candidates can be found in some well-motivated extensions of Standard 
Model. The best-known example would be Minimal Supsersymmetric Standard Model 
(MSSM) with R-parity in which the lightest supersymmetric particle (LSP) turns out to 
be a good DM particle whose thermal relic density is naturally in the right range fl]] . 

Low-energy supersymmetry at the TeV scale has been motivated as a natural expla- 
nation of the electroweak symmetry breaking, and it may also be related to the orgin of 
neutrino masses. This motivates us to introduce into MSSM a TeV-scale seesaw mecha- 
nism explaining the observed neutrino masses and mixing. In this regard, it is interesting 
to address a question whether there can be a new dark matter candidate coming from the 
neutrino sector. 

The most popular way to realize the seesaw mechanism [Q] is to introduce massive 
singlet (right-handed) neutrinos (type I seesaw). In this case, the lightest right-handed 
sneutrino can be stable and become an additional candidate of dark matter whose relic 
density is determined non-thermally || or thermally if the right-handed neutrinos carry 
extra U(l)' interaction [|| or Yukawa interaction associated with the Higgs sector ||. 
Another option is to assume the presence of an SU{2)l triplet with U(l)y charge one 
(Y = 1), which couple to two lepton doublets and thus generate neutrino masses through 
its small vacuum expectation value (type II seesaw). This model predicts quite distinctive 
collider signatures through which the structure of the neutrino mass matrix can be explored 
§. However, its dark matter candidate, the neutral component of the triplet fermion, 
directly couples to the Z gauge boson and thus its large coupling is already ruled out 
by direct detection experiments. The last option of our interest is to introduce Y = 
triplets which couple to lepton and Higgs doublets (type III seesaw) 0. A neutral scalar 
component of such triplet superfields can be the LSP. Let us remark that such a scalar 
triplet LSP can be another realization of the minimal dark matter studied in Ref. §§. 

In the following, we explore properties of this scalar triplet present in the supersym- 
metric type III seesaw mechanism as a dark matter candidate and its signatures in collider 
experiments. 
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2. Triplet spectrum in supersymmetric type III seesaw 



Type III seesaw mechanism introduces real SU(2)l triplets E = (E + , E , £ ) with Y = 0, 
which allows the superpotential, 

Wui = VijLi^Zj + ^MfcSfcEfe . (2.1) 

Integrating out the heavy triplet fields one obtains the seesaw neutrino mass matrix, 

o 2 

M"j = VikVjk^- , (2.2) 



3 



where V2 = {H®). From now on, the generation index for the triplets will be suppressed as 
we are interested in the one containing the lightest neutral scalar component. The fermion 
triplet components, denoted by E^' , have the common mass M. For the mass spectrum 
of the scalar triplet components, denoted by E ±,C) , we need to consider the supersymmetric 
mass M, the soft supersymmetry breaking masses and the electroweak mass splitting. The 
neutral scalar components with T3 = take the masses given by 

m 2 =M 2 + m 2 ± BM (2.3) 

^2,1 

where m is the diagonal soft mass and B is the bilinear soft term which will be assumed to 
be positive without loss of generality. The charged scalar components S ± carrying T3 = ±1 
have the following mass-squared matrix: 



M 2 + fn 2 + c^m 2 z C2{i BM 

BM M 2 + rh 2 — (^m 2 z c 2 ii 



(2.4) 



where c\y is the cosine of the weak mixing angle and the angle (3 is defined by tp = i^/Vl- 
Diagonalizing Eq. ([2.4), one finds the mass-squared eigenvalues, 



77) 



|± =M 2 + rh 2 ± yj B 2 M 2 + c^m%(? 2p . (2.5) 



Note that the lighter scalars, and S^, can have smaller masses than their fermionic 
partners and the ordinary superparticles < M, rh) due to a large negative term 

driven by the .B-term, and can be even lighter than the ordinary lightest supersymmetric 
particle (OLSP) in the MSSM. The tree-level mass splitting between E 1 and driven by 
the electroweak symmetry breaking is given by 

2 4 2 

Am tree = [m £ ± - m t o\ tvce « ~ ± BMm ? ( 2 - 6 ) 

in the limit of m 2 z <C BM. Note that the tree-level mass splitting is negative. On the 
other hand, one-loop correction induces positive contribution given by 

Am loop = -^-±[f( rw ) - c 2 w f(r z )] (2.7) 
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Figure 1: Contour plot of mass splitting Am between the charged and the neutral (DM) scalar 
components as a function of the DM mass ttidm and the parameter b = v BM/wdm- The black 
region is for Am < and thus excluded for our discussion. In the next darkest region, we have 
Am < 100 MeV, and Am is increased by 10 MeV for each contour region. In the brightest region, 
Am is larger than 160 MeV but below the upper limit of 167 MeV. 

where r\y y z = mw^lmff* anc ^ f( r ) = — r [2r 3 lnr + (r 2 — 4) 3 / 2 ln(r 2 — 2 — ryV 2 — 4) /2]/4 
§■ 

The one-loop splitting reaches the typical maximum value Ami oop ~ 167 MeV for 
TUf. S> mw,z which can be partly canceled by the tree- level contribution. In Fig. 1, we 
show the total mass splitting Am = Am loop + Am tree in the plane of the lightest neutral 
scalar mass todm = m^o and the dimensionless parameter b = V BM /m^M quantifying 
the size of the B-term. Apart from the black region in the lower left corner, the neutral 
scalar component remains lighter than the charged scalar and can be the LSP dark matter. 
Depending on the size of Am, some interesting collider signatures will occur as will be 
discussed in Section 4. 



3. Relic density and indirect detection limits 

The dominant annihilation channel of the dark matter, Y,®, is the 'direct' gauge coupling 
with W^, and thus it has a large s-wave annihilation rate [^: 

Ana 2 . 

™ 3.1 

m| ? 

It can give rise to the right thermal relic density only in the multi-TeV range in the standard 
cosmology. For smaller mass, the annihilation cross-section becomes large to suppress the 
standard thermal relic density. In this case, the correct amount of dark matter relics must 
come from a certain non-thermal origin, or from a non-standard thermal history of the 
universe. An interesting possibility for the latter is a kinetic energy (kination) domination 
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Figure 2: The solid line shows the reduced kination-radiation equality temperature, T r = 
\J 9*eql g*f T eq , required to produce the right dark matter density. The dotted line is he freeze- 
out temperature Tf. 



driven by the evolution of quintessence which forms dark energy today ||. Such a kination 
era would immediately follow the period of inflation in a cosmological scenario identifying 



the inflaton and quintessence field [ 10 1 , where one can still find sufficient reheating through 
gravitational particle production |pH| . 

For a phenomenological consideration of kination cosmology, it is enough to identify 
one free parameter T eq at which the kination and radiation components become equal. At 
a given temperature T, the expansion parameter can be expressed as 



H(T) = hJi + ^-^ (3.2) 

V 9*eq 1 eq 



where H st is the Hubble parameter in the standard cosmology, H st ~ l.66^/~g^T 2 /m p i. 
If dark matter freeze-out occurs before T eq (when H > H s t), a larger annihilation cross- 
section is required to produce a right thermal relic density. Given the dark matter freeze- 
out temperature Tf and the kination-radiation equality temperature T eq , one can find the 
approximate formula for the present dark matter density: 

u2 0.88 x 10- 10 GeV- 2 z f z r 
S2dm« ~ ; — ; r=k( — ), (3.3) 

where Zf = m^o/T/, z r = y/g*f / g* eq m^o / T eq , and k(u) = u/hx{u + v 7 ! + u 2 ) @. Note 
that the function k characterizing the kination domination recovers the standard value k = 
1 in the limit of u — > 0. For the purpose of our investigation, we calculate Tf approximately 
by equating the annihilation cross-section and the expansion parameter. Then, requiring 
^dm/i 2 = 0.11, one can find appropriate values of the freeze-out temperature Tf and the 
redefined kination-radiation equality temperature T r = \/g* e q/g*f T eq as functions of the 
dark matter mass mg - The result is shown in Fig. 2. One finds that the required T r is 
larger than about 5 MeV for the dark matter mass larger than 200 GeV, and thus there 
is no conflict with the standard big-bang nucleosynthesis in the parameter region of our 
interest. 



- 4 - 



1 



200 400 600 800 1000 1200 
m DM in GeV 

Figure 3: The cross-section of the scalar triplet dark matter annihilating into W + W~ is shown 
by the thick solid line in unit of 3 x 10 -26 cm 3 /sec. Thin solid and dashed lines summarize var- 
ious astrophysical and cosmological constraints. Solid lines from the uppermost one come from 
the observations of extra-galactic diffuse gamma background [Belikov-Hooper Q], big-bang nucle- 
osynthesis |l6| , anti-proton flux in cosmic rays , and galactic center radio observation (assuming 
NFW profile) |fj|. The upper (lower) short-dashed line is the limit from CMB measurement by 
WMAP5 (PLANCK) The long-dashed line is the optical depth bound [Cirelli et.al. ||]. 



Concerning the direct detection of our dark matter particle, let us recall that it has 
Q = Y = and its leading contribution to the scattering on nuclei comes from scalar 
interaction at one-loop. Thus its nucleonic cross-section is quite small: cr n ucieon ~ 10 -45 cm 2 
Q, which is about two orders of magnitude below the current direct detection limit. On the 
other hand, it is amusing to note that our scalar dark matter triplet with Am ~ 10 MeV 
can explain the DAMA/LIBRA modulation results |l3| through the element-dependent 
resonant scattering 



The dark matter annihilation rate (3.1) of — * W + W~ is much larger than the 
standard freeze-out value: 3 x 10~ 26 cm 3 /sec for sub-TeV dark matter mass, as shown by 
the thick solid line in Fig. 3. Such an enhanced annihilation can lead to various indirect 
signals in astrophysics and cosmology, and thus gets restricted by current observational 
data. Fig. 3 summarizes various constraints analyzed in the recent literature considering 
the effects of dark matter annihilation to the extra-galactic diffuse gamma-ray background 
by produced energetic electrons and positrons |]15| , the changes in light element abundances 
predicted successfully by the big-bang nucleosynthesis |jlq ], the anti-proton flux in cosmic 
rays jl?]], the radio observation from the galactic center fl^| , the extra energy injection 
during the recombination epoch , and the observed optical depth of the Universe p(| . 
The most stringent limit comes from the galactic center radio observation which strongly 
depends on astrophysics models. The NFW profile is assumed for the lowest solid line 
which puts a strong bound: m^o > 900 GeV . But, this bound disappears if the isothermal 
profile is assumed. Therefore, we do not take this bound for our consideration of collider 
signatures. At the moment, the most stringent bound comes from the cosmic ray anti- 
proton flux measurements; 

m g0 > 520 GeV, (3.4) 



-'i 



which is comparable to the future CMB limit. Independently of such indirect detection 
constraints, collider experiments will be able to provide another limit on the scalar triplet 
mass. 



4. LHC signatures 

Some of interesting collider signatures of neutrino mass models at TeV scale can occur 
in lepton-flavor violating decays of new particles responsible for generating light neutrino 
masses. They usually involve displaced vertices related to small neutrino masses and same- 
sign dilepton final states. In type III seesaw, the production of triplet fermions through 



off-shell W boson can yield such a signal: pp -> S ± S° -> l^W^Z f|| || ||, ||. The 
decay rate of the triplet fermion is given by ~ h 2 M/8tt where the neturino Yukawa 
coupling can be traded with the effective neutrino mass: rh v = y 2 v\jM. Thus its lifetime 
becomes rs ~ 0.1 mm for rh u = 0.05 eV and M = vi ~ 174 GeV. This result is applicable 
also to type I seesaw. In type II seesaw, where doubly charged Higgs boson decays to two 
same-sign leptons, the maximum decay length can be found to be about 0.3 mm f6j. 

In our DM scenario in supersymmetric type III seesaw, we concentrate on the collider 
phenomenology of the triplet scalars. A peculiar feature of displaced vertices reflecting 
small neutrino Yukawa couplings occurs now in the decay of the ordinary lightest super- 
symmetric particle. For our discussion, we will take the usual lightest neutralno Xi (mostly 
Bino) as the OLSP. Note that the seesaw mechanism in the fermion (scalar) sector induces 
mixing between lepton (slepton) doublets and fermion (scalar) triplets. In particular, two- 
body decays of the OLSP can arise through the mixing between sleptons and scalar triplets 
induced by the soft A-teim and supersymmetric F-term, 

V mix = m D A [-L + fE+] + m D M [-L ut°* + ZXT*] + h.c. , (4.1) 
V 2 v 2 

where uid denotes the Dirac neutrino mass between v and S°: uid = yv2- The OLSP 
decay rates, for x? — ► J and z^S?, are found to be 

/ 2 \ 2 



16 " Xi \ m 2 

\ xi/ 

r(x!-^) = ^. V (i-^|J ■ (4.2) 

Here the small parameters 0r~ quantify the mixing between (y) and Tif (S?) : 

^ M(M - A) 

Tu ~ 2 2 
1,1/ Ei 




5 where S = \ . (4-3) 



Recall that 5 is the small mixing angle between v (I) and S° (S~) coming from the seesaw 
relation m u = m 2 D /M. Compared to the fermion triplet decay rate Te discussed above, the 
OLSP decay rate ( |4.2|) , r x o ~ a'5 2 m x o/8, is typically smaller by an order of magnitude, 
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and thus we expect to have more distinct displaced vertices. Taking again 6j~ ~ 5 and 
m v = 0.05 eV, the OLSP decay length becomes 



a' _ ™ x oi 1 



M 



0.3 cm (4.4) 



for m^o = M. 

The charged and neutral scalar fields present in our model are nearly degenerate with 
each other. As shown in Section 1, the tree-level mass splitting driven by the B-term and 
D-term masses always gives negative contribution to Am = m^± — m^ , which can partly 
cancel the positive one-loop correction. As a consequence of this, the mass splitting Am 
can be smaller than the maximal splitting of Am ~ 167 MeV induced solely by one-loop 
correction Q. Depending on whether Am is larger or smaller than m n ±, one finds the 
allowed decay channel, whose rates are given by 



T(Et - S?^) = l G^ d Am 3 f^l - ^ ; 

T(tf - t^v) = G 2 F Am 5 . (4.5) 



Note that we have a similar situation as in the Wino LSP scenario [25] with much more 
restricted value of Am. Independently of the scalar triplet mass, the two-body decay rate 
takes its maximum value for Am ~ 167 MeV, and its minumum value for Am = m n ±, 
below which the three body decay channel opens up. Unless the mass splitting Am is finely 
tuned to be less than m e ± = 0.5 MeV, the three body decay rate of Eq. (f4.5|) is applied for 
Am < m w ±. Thus, we restrict ourselves to the decay length of the charged scalar triplet 
in the range: 

5.5 cm < T-i <6.3m. (4.6) 

Recall that ir^ or e ± from the decay is too soft to be observed in the LHC detectors. 
Thus one can only observe highly-ionizing tracks caused by Y,f which typically disappear 
somewhere inside a detector. When such a heavy charged particle moves slowly (/? < 1), 
its momentum and velocity (and thus mass) can be measured simultaneously, providing 
a unique background- free and model- independent search for new physics. Experimental 
studies on the mass reconstruction have been performed by using the time-of-flight mea- 
surement at the muon detectors of ATLAS and CMS, and also by the ionization-rate 
(dE/dx) measurement at the tracker of CMS [^Tj. Both methods can be applied to our 
case if the mass splitting is close to or below m n ± [see Q4.6D ]. 

The heavy scalar triplet can be produced following the (cascade) production of the 
OLSP: 

or directly through W / Z boson exchange: 

PP - { ~ X 1 1 ■ (4.8) 



-7- 



In the first case, one can have both the same-sign dileptons and same-sign anomalous 
charged tracks resulting from the Majorana nature of the mother particle. Note that the 
OLSP mass can be reconstructed from the energy (momentum) measurement of the charged 
leptons and the mass measurement of £ x . Furthermore, the reconstructed displaced ver- 
tices of the OLSP decay (O) can be useful to distinguish the candidate events. The 
production cross-section of the process ( [4.7D depends on models of supersymmetry break- 
ing and sparticle mass spectrum. Assuming the squark and gluino have the same mass, 
their production cross-section at the LHC with yfs = 14 TeV is 2 pb for the squark/gluino 
mass 1 TeV pq] , providing a large number of the lepton number violating processes and 
anomalous charged tracks ( |4.7| ) even at the early stage of the LHC. From the result of 
Ref. J28|, we can infer that the scalar triplet signals with a few events can be probed up to 
the squark/gluino mass of 1.7 (2.8) TeV for the luminosity of 0.1 (10) fb . 

Apart from this, the process ( ^4.8|) provides another means to test our model, inde- 
pendently of the squark and gluino mass. The corresponding cross-sections have been 
calculated in Ref. [26]. For = 550 GeV, for instance, the production cross-sections 
for S^S? and are 2.2 fb and 1.1 fb, respectively. Thus, copious (single and double) 

anomalous charged tracks can be observed if the dark matter mass is close to its lower bound 
and the integrated luminosity reaches a nominal value of 10 fb -1 . The largest scalar triplet 
mass, for which the scenario under discussion can be probed at the LHC luminosity of 10 
(100) fb -1 , is about 800 (1250) GeV corresponding to the total cross-section of 0.5 (0.05) 
fb. We note that, given the dark matter mass bound of 520 GeV, the direct production 
is too small to be observed for the initial luminosity of 100 pb -1 . However, the cascade 
production from the initial squark/gluino decay could be observed if the squark/gluino 
mass is smaller than about 1.7 TeV as mentioned above. 

Let us finally make a remark on the fermion triplet signature. When the fermion mass 
M is not too large, one can look for the same-sign dilepton signal pp — > l^l^ + 4j followed 
by the E ± S° production as mentioned earlier in this section. Note that its production rate 
is 10 times larger than the scalar production rate, and the above lepton number violating 
signals can be probed up to a fermion mass 450 (700) GeV at the LHC luminosity of 
10 (100) fb -1 [p3|. Let us also recall that the trilepton signal Z ± Z ± Z =F is as good as the 



same-sign dilepton signal in probing the fermion triplet [22]. 



5. Conclusion 

It is suggested that the neutral scalar component °f the Y= triplet in the supersym- 
metric type III seesaw mechanism can be a viable dark matter candidate. The 5-term 
contribution can lead to a large mass splitting between two neutral scalars and the lighter 
one becomes the LSP. For its mass in the sub- TeV range, the annihilation cross-section is 
much larger than the standard value 3 x 10~ 26 cm 3 /sec required for a thermal production 
of the observed relic density. Such an enhanced annihilation can be consistently accom- 
modated by considering a non-standard cosmology: kination domination around the dark 
matter freeze-out era. This kind of scenario could be tested or limited indirectly through 
various astrophysical observations as summarized in Fig. 3. Currently the strongest bound 
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comes from the radio observation from the galactic center which is however strongly de- 
pendent on the dark matter profile. Next stringent constraint is put by the cosmic ray 
anti-proton measurement, requiring m^o > 520 GeV. 

A peculiar aspect of our model in view of collider phenomenology is the existence of 
the charged scalar particle Sf which is almost degenerate with the dark matter particle Y,®. 
Their mass splitting Am gets both the (negative) tree-level and (positive) one-loop con- 
tributions which can cancel each other. When Am < m w ±, the charged scalar allows only 
three-body decay, — > Y^e^v, with r ~ 6.3 m. This charged partner of the dark mat- 
ter is another example of heavy long-lived charged particles whose mass can be measured 
model-independently by detecting their anomalous tracks in the LHC experiments. The 
charged scalars can be produced directly through the W/Z boson exchange, or after the 
production and the decay of the ordinary lightest supersymmetric particle like the Bino. 
In the latter case, we can probe clean signatures of same-sign dilepton and di-charged- 
scalar tracks having measurable displaced vertices (r ~ 0.3 cm). Furthermore, the mass 
measurement of the charged scalars allows us to determine the OLSP mass as well. 

Given the direct production cross-section of our scalar particles, pp — > S^S? and 
TjfYi, numerous anomalous charged particle tracks could be observed if the scalar dark 
matter mass is in its low end; e.g., about 30 events for the mass 550 GeV at the luminosity 
10 fb _1 . The discovery /exclusion of the model can be made up to m^o < 1250 GeV 
with the integrated luminosity 100 fb _1 as no background events are expected. On the 
other hand, the previous studies showed that the direct production of the fermion triplet 
pp — > leading to the same-sign dilepton signals can be probed up to the fermion 

triplet mass 700 GeV. Of course, one could observe much more events if a squark or gluino 
and thus the OLSP can be copiously produced, which depends on the mass spectrum of 
heavier supersymmetric particles. 

Acknowlegement: The author thanks Suyong Choi and Manuel Drees for valuable 
comments on the LHC signals. This work was supported by Korea Neutrino Research 
Center through National Research Foundation of Korea Grant (2009-0083526). 

References 

[1] For a review, see, G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 
(1996) [arXiv:hep-ph/9506380]. 

[2] For a review, see, R. N. Mohapatra et al, "Theory of neutrinos: A white paper," Rept. Prog. 
Phys. 70, 1757 (2007) [arXiv:hep-ph/0510213]. 

[3] T. Asaka, K. Ishiwata and T. Moroi, Phys. Rev. D 73, 051301 (2006) [arXiv:hep-ph/0512118]. 

[4] H.-S. Lee, K. T. Matchev and S. Nasri, Phys. Rev. D 76, 041302 (2007) 
[arXiv:hep-ph/0702223]. 

[5] D. G. Cerdeno, C. Munoz and O. Seto, Phys. Rev. D 79, 023510 (2009) [arXiv:0807.3029 
[hep-ph]]; F. Dcppisch and A. Pilaftsis, JHEP 0810, 080 (2008) [arXiv:0808.0490 [hcp-ph]]. 

[6] E. J. Chun, K. Y. Lee and S. C. Park, Phys. Lett. B 566, 142 (2003) [arXiv:hep-ph/0304069]. 

[7] R. Foot, H. Lew, X.-G. He and G. C. Joshi, Z. Phys. C44 (1989) 441. 



- 9 - 



[8] M. Cirelli, N. Fornengo and A. Strumia, Nucl. Phys. B 753, 178 (2006) 
[arXiv:hep-ph/0512090]. 

[9] R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80 1582 (1998) 

[arXiv:astro-ph/9708069]; P. J. Steinhardt, L. M. Wang and I. Zlatev, Phys. Rev. D 59 
123504 (1999) [arXiv:astro-ph/9812313]. 

101 B. Spokoiny, Phys. Lett. B 315 40 [arXiv:gr-qc/9306008]; P. J. E. Peebles and A. Vilenkin, 
Phys. Rev. D 59 063505 (1999) [arXiv:astro-ph/9810509]. 

I ll L. H. Ford, Phys. Rev. D 35, 2955 (1987); E. J. Chun, S. Scopel and I. Zaballa, JCAP 0907, 
022 (2009) [arXiv:0904.0675 [hcp-ph]]. 

121 P. Salati, Phys. Lett. B 571, 121 (2003) [arXiv:astro-ph/0207396]. 

131 Y. Bai and P. J. Fox, arXiv:0909.2900 [hep-ph]. 

I tl M. Pospelov and A. Ritz, Phys. Rev. D 78, 055003 (2008) [arXiv:0803.2251 [hcp-ph]]. 

I 51 M. Kawasaki, K. Kohri and K. Nakayama, Phys. Rev. D 80, 023517 (2009) [arXiv:0904.3626 
[astro-ph.CO]]; S. Profumo and T. E. Jeltema, JCAP 0907, 020 (2009) [arXiv:0906.0001 
[astro-ph.CO]]; A. V. Belikov and D. Hooper, arXiv:0906.2251 [astro-ph.CO]. 

161 J- Hisano, M. Kawasaki, K. Kohri, T. Moroi and K. Nakayama, Phys. Rev. D 79, 083522 
(2009) [arXiv:0901.3582 [hcp-ph]]. 

171 F. Donato, D. Maurin, P. Brun, T. Delahaye and P. Salati, Phys. Rev. Lett. 102, 071301 
(2009) [arXiv:0810.5292 [astro-ph]]. 

LSI G. Bertone, M. Cirelli, A. Strumia and M. Taoso, JCAP 0903, 009 (2009) [arXiv:081 1.3744 
[astro-ph]]. 

1.91 S. Galli, F. Iocco, G. Bertone and A. Melchiorri, Phys. Rev. D 80, 023505 (2009) 

[arXiv:0905.0003 [astro-ph.CO]]; T. R. Slatyer, N. Padmanabhan and D. P. Finkbcincr, 
arXiv:0906.1197 [astro-ph.CO]. 

201 G. Huetsi, A. Hektor and M. Raidal, arXiv:0906.4550 [astro-ph.CO]; M. Cirelli, F. Iocco and 
P. Panci, arXiv:0907.0719 [astro-ph.CO]. 

211 R. Franceschini, T. Hambye and A. Strumia, Phys. Rev. D 78, 033002 (2008) 
[arXiv:0805.1613 [hcp-ph]]. 

221 F. del Aguila and J. A. Aguilar-Saavedra, Nucl. Phys. B 813, 22 (2009) [arXiv:0808.2468 
[hcp-ph]]. 

231 A. Arhrib, B. Bajc, D. K. Ghosh, T. Han, G. Y. Huang, I. Puljak and G. Senjanovic, 
arXiv:0904.2390 [hep-ph]. 

24] T. Li and X. G. He, arXiv:0907.4193 [hep-ph]. 

25] C. H. Chen, M. Drees and J. F. Gunion, Phys. Rev. Lett. 76, 2002 (1996) 

[arXiv:hcp-ph/9512230]; J. L. Feng, T. Moroi, L. Randall, M. Strassler and S. f. Su, Phys. 
Rev. Lett. 83, 1731 (1999) [arXiv:hep-ph/9904250]. 

[26] P. Fileviez Perez, H. H. Patel, M. J. Ramsey-Musolf and K. Wang, Phys. Rev. D 79, 055024 
(2009) [arXiv:081 1.3957 [hep-ph]]. 



- 10 - 



[27] S. Giagu, "Search for long lived particles in ATLAS and CMS", Talk presented at ICHEP 
2008, Philadelphia, PA; The ATLAS Collaboration, "Expected Performance of the ATLAS 
Experiment : Detector, Trigger and Physics", CERN-OPEN-2008-020 [arXiv:0901.0512 
[hep-ex]]; The CMS Collaboration, "Search for heavy stable charged particles with 100/pb 
and 1/fb in the CMS experiment", CMS-PAS-EXO-08-003. 

[28] W. Beenakker, S. Brcnsing, M. Kramer, A. Kulesza, E. Laenen and I. Niessen, 
arXiv:0909.4418 [hep-ph]. 



- 11 - 



